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The "canonical" variables of the Kosterlitz-Thouless theory — fields $o(r) and fir), generally 
believed to stand for vortices and phonons (or their XY equivalents, like spin waves, etc.) turn 
out to be neither vortices and phonons, nor, strictly speaking, canonical variables. The latter fact 
explains paradoxes of (i) absence of interaction between $o and if, and (ii) non-physical contribution 
of small vortex pairs to long-range phase correlations. We resolve the paradoxes by explicitly relating 
"l>o and LP to canonical vortex-pair and phonon variables. 
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Three decades ago, Kosterlitz and Thouless developed 
an accurate renormalization-group description of what is 
now called a Berezinskii-Kosterlitz-Thouless (BKT) tran- 
sition 0, Q — a phase transitions in a wide class of two- 
dimensional systems characterized by short-range inter- 
actions and global U(l) symmetry. In accordance with 
the Mermin- Wagner theorem Q, such systems can not 
exhibit long-range order at any finite temperature. In- 
stead, the low-temperature phase features divergent long- 
wave fluctuations resulting in a power-law decay of phase 
correlations at large distances Kosterlitz and Thou- 
less revealed the importance of configurations with point- 
like topological defects or topological charges, such as 
Coulomb charges in plasma, dislocations in crystals, vor- 
tices in superfluid and spin systems, etc. At low tem- 
perature, the defects exist only in the form a dilute gas 
of bound pairs of opposite topological charges. At higher 
temperature, pairs of large separation become more prob- 
able and eventual pair dissociation at critical tempera- 
ture destroys the algebraic long-range order. Extensive 
theoretical work 0,13, IE IS^Q provides a complete 
quantitative description of critical properties in such sys- 
tems. The theory is corroborated by comprehensive ex- 
perimental studies of ''He films H, l^, hll and super- 
conducting Josephson arrays ^2 ■ Recent advances in the 
area of ultra-cold gases have made it possible to render 
BKT transition in optical lattices |3 . 

The Kosterlitz-Thouless (KT) theory starts with a 
generic effective action 



(1) 



For definiteness, we consider the case of a superfluid film, 
in which the field <I>(r) has the meaning of the velocity po- 
tential, V — (l/m)V<I' (we set h = 1), and Kq = no/2mT, 
where m is the mass of a particle and uq is the "bare" 
superfluid number density obtained by averaging out mi- 
croscopic fluctuations up to some mesoscopic scale Iq. 
Hence, no = "-oCo)- The field (f> is then split, 

$ = $0 + , (2) 

into a singular part $0j containing all the topological 
defects, and a regular part (p. 



By definition, <i>o satisfies the non-zero velocity circu- 
lation condition 



f V$o(r) 



■ dr = 2nlj, 



(3) 



where Cj is a contour enclosing only the j-th defect and 
Ij are integers, while V(/5 is circulation-free. The next 
crucial step is to require 

A$o(r) - (4) 

(except for the isolated points of defects). The standard 
motivation of Eq. Q is that it guarantees that $0 mini- 
mizes the action when = Q. The definitions of $0 and 
(f thus become unambiguous and, most importantly, the 
action takes the form of two independent terms: 



(5) 



At this point, one conjectures that $0 f^nd ip corre- 
spond to vortices and phonons (spin waves, etc.), respec- 
tively. This identification, which might seem to be quite 
natural — or at least merely terminological and mathe- 
matically irrelevant — is not that innocent. If the two 
fields are not canonical vortices and phonons, then one 
faces a problem of justifying writing the partition func- 
tion in the form 



oc 



/ exp{-A[^]}I?^ / exp{-A[$o]}n^^ 



(6) 



where rj is the position of the j-th defect and N is the 
number of defects. This expression should also include 
the Jacobian of the transformation from canonical vari- 
ables. Remarkable agreement between the Kosterlitz- 
Thouless theory and experimental data suggests that the 
Jacobian is unimportant, but without explicitly demon- 
strating this fact the theory is incomplete. 

Apart from Jacobian, there is also an issue of peculiar 
"collective" behavior of formally independent fields $0 
and (f. Asymptotic long-range phase correlations in a 
superfluid are due to phonons. The corresponding action 
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in terms of the genuine phonon field (p is 



2mT 



(7) 



with Us the macroscopic superfluid density. Would (p 
stand for phonons, we were to identify its long-wave har- 
monics with ip. This, however, would imply tiqCo) — "-s, 
which is definitely not the case since $0 7^ at the length 
scale Iq. This paradox can be formulated as an observa- 
tion that it is impossible to renormalize the sound veloc- 
ity by vortex pairs without the phonon-vortex coupling. 
The only logical solution is then that ip is not a phonon 
field. 

Another paradoxical circumstance is associated with 
interpreting $0 a purely vortex field. In a 2D super- 
fluid, all vortices are bound in microscopic pairs and one 
would not expect them to be directly observable in long- 
range correlation properties. The only physical way for 
vortex pairs to manifest themselves at the macroscopic 
scale is to renormalize the superfluid density. However, 
the requirement 10} implies that vortex pairs do con- 
tribute to the long-range correlations. The way they do 
it reveals a conspiracy between $0 and (p Q. Due to 
statistical independence of the fields $0 and ip, the one- 
particle density matrix at large distances, 

p(r) cx ( exp[i$(r) - i$(0)] ) , (8) 

factorizes: p{r) oc r$f, (r) r;^(r), where 

r*„(r) = (exp[^$o(r)-^$o(0)] ) , (9) 



r^(r) = ( expMr) - z^(0)] ) 



(10) 



Remarkably, the independent correlation functions F^^ 
and make no physical sense separately, since they both 
depend on the bare superfluid density uq, but when the 
two are combined in p(r), the parameter uq drops out, 
and the density matrix decays with the proper exponent 
mT/2nns 4]. 

A rational explanation of this "secret agreement" be- 
tween $0 and ip is that being mathematically indepen- 
dent, the two fields are deeply connected physically. The 
above-mentioned structure of the correlation function 
even suggests a qualitative form of the connection: The 
long- wave part of the vortex-pair field $0 actually belongs 
to phonons, not vortices. 

In what follows, we trace the model ^ back to its 
dynamical Hamiltonian form and derive the canonical 
parametrization of vortices and phonons from that start- 
ing point. In doing so, we utilize the formalism recently 
developed by two of us 0|, from which it is directly 
seen that the positions of vortices {rj} and the field (p 
are canonical variables only in the limiting cases of in- 
compressible fluid {ip = 0) and vortex-free environment 
($0 = 0) respectively. An explicit transform from {fj} 
and ip{r) to the canonical variables justifies that the de- 
viation of the Jacobian from unity is irrelevant in the 



context of the KT theory. Our analysis allows us to refor- 
mulate KT theory in terms of canonical variables: vortex 
pairs and genuine phonons, coupled to each other in the 
most intuitive way: Vortex pairs interact with the long- 
wave fluctuations of the velocity field precisely the same 
way they interact with a homogeneous velocity field. This 
interaction naturally accounts for the renormalization of 
the long-wavelength-phonon stiffness and leads to the 
coarse-grained effective action in the form of Eq. 0. 
In complete agreement with physical understanding of 
vortex pairs as essentially local objects, we demonstrate 
that the far- field of ^o{r) belongs to phonons. Corre- 
spondingly, the long-range decay of phase correlations 
is governed by the statistics of long-wavelength phonons 
only, i.e. by the effective action ||7J). 

Putting aside the issue of the Jacobian and pior to the 
discussion of the dynamic model, a purely statistical in- 
sight into the problem can be obtained by constructing 
an alternative to {^o^ip) set of variables. [For simplic- 
ity, below we deal with only one vortex-anti-vortex pair; 
the generalization to finite (but small) density of pairs 
is straightforward.] Consider a vortex pair of separation 
R = ri — r2 located at the point rp = (ri -|- r2)/2, where 
ri and r2 are the positions of the vortices with li = 1 
and h = respectively. Introduce an auxiliary field 
(po{r) such that it approaches $o(r) when |r — rp| » R, 
and, in contrast to <I'o(r), is regular at all distances. This 
definition fixes the dipole moment of the new field: 



rA(^o(r) 



27r(R X z) , 



(11) 



where z is a unit vector along z-axis. Make a transfor- 
mation 



(p{r) = p(r) - (po{r) , 



(12) 



which just shifts the field by a regular (rp,R)- 
depcndent field ipo, and thus does not change the con- 
figurational volume: "Dip = Dtp. After this transforma- 
tion, the long-range behavior of the density matrix is 
completely described in terms of the filed (p: 



p{r) oc ( exp[i.^(r) - i<p{0)] ) 



(13) 



This simplification comes at a price: the vortex pair now 
couples to (p. The structure of the interaction term be- 
tween the pair and the long-wave harmonics of the field 
(p (such that X ^ R, where A is the characteristic wave- 
length) is most transparent. It reads 



A, 



27rno 
mT 



(Rx z) ■ Vip 



(14) 



i.e. the vortex pair interacts with the long-wave part of 
ip exactly the same way it would interact with a homo- 



geneous velocity flow (1/m) 



One does not have 



to take this interaction into account explicitly in the KT 
renormalization group treatment, since its only relevant 
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effect is to replace Uq with Ug for phonons. [The effect 
of phonons on the statistics of vortex pairs is neghgibly 
small in the limit of i? ^ oo, as is clear from a direct 
estimate, see also below.] There is little doubt at this 
point that (p corresponds to genuine phonons and one 
just needs to formally demonstrate this fact. 

We start with Popov's hydrodynamic Lagrangian 



L 



H= d 



H 



"0 Iy7 |2 



1 

2m: 



V 



+ ?-|V$o|^ + "ovo • V4>o 



(15) 



(16) 



Here, the energy functional H has been expanded to the 
leading order with respect to small density fluctuations 
r) <C no, $0 and (p are defined by Eqs. Q-Q, vq is the 
velocity of a global flow, and x is the compressibility. 
The typical vortex core size ~ ao = y^>c/nQm is much 
smaller than any other physical length scale. 

If the term J (Pr tj^q = T were absent, rj and (p would 
be the canonical conjugate phonon variables, while 



d rrio$o 



-2nno IjUji. 
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3 ; 



(17) 



where (xj 



y3i 



Tj, would imply that Xj and j/j are the 



canonical conjugate vortex variables. However, T is lin- 
ear in the time derivatives of Xj and yj and also contains 
r] making the set of variables {77, (p}, {vj} not canonical, 
and meaning that H in these variables is not a Hamilto- 
nian. 

We are interested here only in the KT theory for the su- 
perfluid phase in the vicinity of the transition, including 
the critical point, where the concentration of vortex pairs 
of size ~ i? is much smaller than as i? — > 00. Cor- 
respondingly, at a phonon wavelength A only pairs with 
R X contribute to the renormalization of the sound 
velocity. This allows us to use the small parameter 



-R/A< 1 



(18) 



for deriving canonical variables in the form of a regular 
perturbative expansion starting from the zeroth approx- 
imation {77, (/s, r^} 14]. 



It is straightforward to show that 



r = ^2^?, [zx VQ(r,)].r, , 



(19) 



where Q{r) is defined by AQ(r) = ri{r). We first switch 
to the Fourier representation of {ry, (p}: 



p{v) = ^ Jll2Vu:^^ [ e^^^ Cq - e'^^^ c* ] , (20) 



where = {yj rioj ycni) q and {cq,c*} are resembling 
(and to the zeroth approximation coincide with) the 
classical-field counterparts of phonon creation and an- 
nihilation operators, and V is the system volume. Let 
the vortex (Ji = 1) and the anti- vortex (Z2 = — 1) in a 
pair have coordinates ri = rp + R/2 and r2 = Tp — R/2 
respectively. Next, we expand QijCj) in T into series with 
respect to qR <C 1. The resulting terms are eliminated 
by iteratively correcting the variables {r^}, {cq,c*} as 
described in Ref. jl3], so that the Lagrangian takes on 
the canonical form 



L = > i c„c* 



q 



L-q<-q 



27r7io I 



3V 3^ 3 



H{r 



3 ' "-qi 



(21) 



where = (xj,yj) and Cq,c* are the Hamiltonian vari- 
ables. For our purposes we shall need only the leading 
correction, which does not change the vortex variables, 

= f, , (22) 

and for the phonon variables, yields 



27r 



m apq [qxRy - qyRx 



o^q""? 



2V2 



q- 



(23) 



where an equivalent set of vortex variables, R = f 1 — f2, 
fp = (fi-|-f2)/2, is used. The Jacobian of the transforma- 
tion (|^ - (|^ equals unity. If higher-order (in 77/no <C 1 
and qR ^ 1) terms are included in Eqs. if^ - lf^ . the 
deviation of the Jacobian from unity is of the order of 
{ri/n()){qR). From now on we omit the tildes over the 
vortex canonical variables in view of Eq. H22|l . 

The canonical phonon fields, 77, (p, are defined anal- 
ogously to if^ . with {cq, ci} replacing {cq, c*}. After 
substituting H23|) for Cq in H2U|) and taking the sum over 
q the original variable p is expressed in terms of the 
canonical variables as 



p{r) = ^(r) - 



[(r - rp)yR^ - (r - rp)^Ry] 



(24) 



Now we note that at large distances from the vortex pair 
r - Tp] > i?, the field $o(r) is given by 0] 



*o(r) 



[(r - rp)yRa, - (r ~ rp)^Ry] 
]r - r„|2 



(25) 



Along with Eq. H24|) , this implies that sufRciently far from 
Tp, $o(r) does not belong to the vortex- anti- vortex pair 
at all, but is actually a part of the phonon field (p. 

After a standard algebra, the Hamiltonian assumes the 
form: 



27rno 



ph 



'd' 



ln(i?//o) + 2E,, 
no I ,2 1 ^2 



(26) 



i?i„ti = 27rno (R X z) • vo, 
27rno 

Hint2 = (R X z) • Vip 
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This form is almost identical to the original effective ac- 
tion in terms of R and of Refs. |3jlin3i0- Besides the 
presence of rj, which is trivially integrated out in the par- 
tition function, the only distinctive feature of the Hamil- 
tonian H26|l is the term -ffint2, which couples the vortex 
dipole moment R to the fluid velocity in the sound wave 
oc Vi^l . 

Consider the thermodynamics of the system H26|l near 
the critical point Tj, = ■KUsjlra. The coupling term -ffint2 
does not change the statistics of vortices, since its typical 
value is small, -ffint2/T' ~ qR ^ 1, whereas the contribu- 
tion of i?v diverges logarithmically. Therefore, the su- 
perfluid density, Hs — (l/mV) d^F/dvQ^\^^_^, a = x,y, 
where F = — TlnZ, is given by the Kosterhtz renormal- 
ization group flow In contrast, i?int2 is essential for 
the phonon statistics. Since its structure is identical to 
the vortex coupling to the uniform flow, averaging over R 
and Tp straightforwardly leads to the coarse-grained ef- 
fective action for the long-wavelength phonons governed 
by the renormalized stiffness, Eq. 0. 

To summarize, we have shown that the simplicity of the 



parametrization 101-© — the statistical independence of 
the fields $0 and ip — comes at a price of substantial lack 
of its physical meaning, apart from inconvenience of cal- 
culating off-diagonal correlators, where direct contribu- 
tion of vortex pairs has to be explicitly evaluated with 
the only goal to replace bare superfluid density uq with 
its renormalized value. An alternative parametrization 
in terms of phonon variables (or their XY equivalents) 
renders the Kosterlitz-Thouless scheme even more math- 
ematically simple and accurate, while making it physi- 
cally transparent. The vortex-phonon interaction that 
appears in the effective Hamiltonian does not lead to any 
complications, because the structure of this interaction is 
exactly the same as that of the interaction of vortex pairs 
with a homogeneous external flow and the only effect of 
this interaction is to ensure that both phonons and vor- 
tices are controlled by the renormalized superfluid den- 
sity. 
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